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Abstract 

Longitudinal ion acceleration from high-intensity (I ~ 10 20 Wcm -2 ) laser interactions with 
helium gas jet targets (n e ~ 0.04n c ) have been observed. The ion beam has a maximum energy for 
He 2+ of (40t 3 ) MeV and was directional along the laser propagation path, with the highest energy 
ions being collimated to a cone of less than 10°. 2D particle-in-cell simulations have been used 
to investigate the acceleration mechanism. The time varying magnetic field associated with the 
fast electron current provides a contribution to the accelerating electric field as well as providing 
a collimating field for the ions. A strong correlation between the plasma density and the ion 
acceleration was found. A short plasma scale-length at the vacuum interface was observed to be 
beneficial for the maximum ion energies, but the collimation appears to be improved with longer 
scale-lengths due to enhanced magnetic fields in the ramp acceleration region. 



I. INTRODUCTION 



Laser intensities currently available can accelerate electrons directly with the laser field 
to energies of many times their rest mass, m e = 511 keV, making them highly relativistic. 
The comparatively heavy ions only become relativistic as a direct consequence of the laser 
field when the normalised vector potential, a , exceeds 1836, which corresponds to IX 2 = 
4.6 x 10 24 Wcm _2 /xm 2 . Therefore experimental investigation is currently confined to non- 
relativistic intensities for ions and their movement is governed by the bulk movement of 
electrons. However, extremely large accelerating gradients (> 10 12 Vm" 1 ) can be supported 
in plasmas, which means that the scale of such plasma based accelerators is reduced by orders 
of magnitude compared with conventional designs and so ions can be readily accelerated to 
non-relativist ic energies with current laser systems. 

The non-relativistic critical plasma density for a laser with a frequency ujl is n c = 
6om e uj1/e 2 and is the plasma density above which the laser is unable to propagate. Above 
this density the plasma is known as overdense and below is underdense. In high-intensity 
laser interactions with solid targets (usually overdense), electrons can be accelerated through 
a number of mechanisms, including vacuum heating [1], j x B heating [5] and direct laser 
acceleration [3j. Relativistic electrons can be generated in this way and as they attempt to 
leave the target and move into the surrounding vacuum a large space-charge electric field 
is generated. Ions can be accelerated to 10s of MeV in a direction normal to the target 
surface. Therefore this is known as the target normal sheath acceleration (TNSA) mech- 
anism. Solid targets are used in the majority of proton and ion acceleration experiments, 
taking advantage of the thin (~ nm) hydrocarbon layers formed on the surface of targets as 
the source of protons [U [5]. The production of high-flux, small-emittance [6] and perhaps 
even quasi-mono-energetic [71 [8] proton beams from laser-plasma interactions have numer- 
ous and diverse potential applications; less invasive radiotherapy [9], isotope production for 
positron emission tomography (PET) [TO] , proton radiography of electromagnetic fields [TT] . 
high brightness neutron sources or to deliver the ignition energy in the fast ignitor inertial 
confinement fusion concept [T2] . 

Underdense plasmas have provided an efficient medium for the acceleration of electrons 
to high energies: direct laser acceleration has produced electrons with energies of up to 
350 MeV with temperatures exceeding the ponderomotive potential [3] and laser wakefield 



2 



acceleration has demonstrated that monoenergetic electron bunches of ^ 1 GeV can be 
generated [13] . Ion acceleration from these underdense plasmas has previously been limited 
to the transverse direction with ion acceleration via the ponderomotive Coulomb explosion 
[T4] and shock acceleration [15] mechanisms. More recently a beam of multi-MeV ions has 
been measured from an underdense plasma with an ultraintense laser pulse [T6l [T7] and this 
paper provides further insight into the interaction and ion acceleration mechanism. The 
experimental results presented in reference [16] are reviewed and followed by particle-in- 
cell (PIC) simulation results investigating the acceleration mechanism and the effect of the 
plasma density and scale-length on the ion acceleration. 

II. EXPERIMENT 

A. Experimental set up 

The experiment was performed using the Vulcan Petawatt laser [18] at the Rutherford Ap- 
pleton Laboratory. The laser pulse had a duration of tl ~ 1.0 ps full- width-half-maximum 
(FWHM) with energy of up to 340 J on target to give a maximum power of 340 TW. The 
central wavelength of the Vulcan Petawatt laser is A = 1.054 /xm so the non-relativistic 
critical density is n c — 1.0 x 10 21 cm -3 . An f/3 off-axis parabolic mirror focused the laser 
to spot with a FWHM diameter of 7 /xm to provide a cycle averaged peak vacuum intensity 
of 1.5 x 10 20 Wcm -2 , which corresponds to a peak ag ~ 21. The contrast ratio was ~ 10 -5 . 

The laser was focused to the edge of the gas flow from a 2 mm diameter supersonic 
nozzle. Supersonic nozzles were used to ensure a uniform density profile over a given distance 
and a reasonably sharp density ramp (< 250 /am) both at the front and back of the gas 
jet. The density profile of the gas jet was determined prior to the high-intensity shots by 
interferometry. The backing pressure could be varied so that the electron density of the fully 
ionised plasma could be set between (0.7-4.0) x 10 19 cm -3 . The frequency separation of the 
forward Raman scattered laser spectra (Au = uo pe = ^/n e e 2 /m e e ) confirmed the electron 
density n e of the interaction. 

Ion energy spectra were taken at four different angles from the laser axis to measure 
the angular emission. Thomson ion spectrometers were placed at 90° and 45° to the laser 
propagation direction. In addition, there were two charged particle spectrometers at 10° and 
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0° to the laser propagation direction. Though primarily to measure the electron spectra, 
the open geometry of the magnetic spectrometers also allowed the measurement of ions and 
other positively charged particles. Ions are deflected in the opposite direction from electrons 
in the magnetic field. The deflection of ions in a magnetic field depends simply on their 
mass-to-charge ratio and their momentum, thus allowing their energy to be determined. The 
nuclear track detector CR39, which is insensitive to positrons, was used to detect the ions. 

For a helium gas target, there are two possible ion species, He 1+ and He 2+ . In the 
Thomson spectrometer, these different charge-to-mass ratio species are separated spatially 
by an electric field in addition to the magnetic field. In the magnetic spectrometers, the 
different species are not separated; therefore, it was necessary to distinguish between the 
species by the pit size for all ions deflected by a similar amount. At a particular point on 
the detector, both ion species may be present but the He 2+ will have more energy. The 
higher the energy of an ion, the further into the CR39 it will travel before it is stopped 
and therefore the damage will be deeper into the material. The CR39 is etched in a NaOH 
solution so that pits are formed in the damaged regions. Smaller pits on the surface are the 
result of deeper damage. Hence at a particular point, using the electron spectrometers, the 
smaller pits will be from the He 2+ and the larger ones from He 1+ ions, allowing the species 
to be differentiated. Electron spectra were measured simultaneously during the shot using 
image plate detectors as reported previously [3]. 

B. Experimental results 

For plasma densities of < 2 x 10 19 cm -3 , energetic ions were found to be emitted primarily 
in the transverse direction consistent with previous measurements [IUCE]. However, above 
this density, a clear, reproducible signal of ions was observed in the 0° and 10° spectrometers. 
Presented here is a shot which had an on-target energy of 340 J (peak I « 6 x 10 20 Wcm" 2 ) 
and was incident on a helium plasma with n e ~ 4 x 10 19 cm -3 . 

On this shot, the on- axis (0°) electron spectrum extended to 66 ± 2 MeV, with a char- 
acteristic temperature of T e « 7.4 MeV. The 10° electron spectrum was a little hotter, 
T e « 9.8 MeV, possibly due to filamentation and hosing in the interaction. In the transverse 
direction, the He 2+ had a maximum energy of (7.8 ± 0.6) MeV and the He 1+ ions had a 
maximum energy of (3.4 ± 0.3) MeV. In the longitudinal direction, the He 2+ had a maxi- 
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FIG. 1: A polar plot showing the number of He 2+ ions at different diagnostic angles for 2 MeV 
and 5 MeV ions. 



mum energy of (40lg) MeV and the He 1+ ions had a maximum energy of (10^) MeV. Fig. 
[j] shows a polar plot of the number of He 2+ ions with 2 MeV and 5 MeV at the different 
diagnostic angles. The He 2+ ions have formed two lobes, one in the forward direction along 
the laser axis and the other in the radial direction. Hence, it is deduced that there is a well 
collimated ion beam in the forward direction with a divergence angle of less that 10°. 



III. SIMULATIONS 



The 2D3V particle-in-cell code Osiris [19] was used to investigate the ion acceleration 
mechanism and the effect of various parameters. The simulations were performed on the 48 
node 'Caesar' cluster and on the CXI supercomputer at Imperial College. Stationary simu- 
lation boxes were used to observe the plasma even after the laser has passed. The stationary 
simulation box was 251 x 251 /im, with a resolution of 20.9 cells /A in the longitudinal (x) 
direction and 12.6 cells /A in the transverse direction (y). Since ionisation effects are not 
included in this version of the code, a He 2+ plasma was used to avoid complications due to 
balance of the different charge states. The ion acceleration was investigated with various 
values of plasma electron density, of 0.01 n C) 0.05 n c , 0.1 n C) 0.3 n C) 0.5 n C) 0.75 n c and 0.9 n c 
where n c is the non-relativistic critical density. The density profile for these simulations had 
a 71 /xm long linear density ramp at the front of the plasma, 25 /im of plasma at maximum 
density and a 71 fim density ramp at the back of the target. The vacuum then extended 
a further 84 /xm behind the plasma. The laser pulse was linearly polarised with the laser's 
electric field in the ^/-direction, with a full-width-half-maximum pulse length of t\ = 500 fs 
and a wavelength A = 1.053 /xm. It was focused to a full- width-half-maximum diameter 
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spot wo = 8 jim at the top of the front density ramp, to give a peak normalised vector 
potential 15. 

Another parameter investigated was the length of the density ramp at the rear of the 
target. In solid target ion acceleration, density ramps formed by pre-plasma at the rear of 
the target have been found to be detrimental to ion acceleration [20]. The reason for this is, if 
the pre-plasma has a scale length of greater than the Debye length, A^ = y 'eofc^T ~j / e 2 n e0 , the 
electrons are able to move to shield ions from the accelerating electric field preventing efficient 
acceleration. For solid targets A^ will be of the order of nanometers, but for the underdense 
targets discussed here, A^ is of the order of micrometers. Taking the electron temperature 
to be T e = 7 MeV and the electron density to be 4 x 10 19 cm -3 gives a ~ 3 /im. It 
is therefore expected that the density ramp on the gas jet targets used in the experiment 
will have a detrimental effect on the ion acceleration, and that the simulations can provide 
an approximate scaling with the ramp length. For these simulations, the simulation box 
was 598 x 168 /im with a resolution of 19.8 cells /A in the longitudinal (x) direction and 
9.4 cells /A in the transverse direction (y). Again, the ion species simulated was He 2+ and 
the initial plasma electron density used was 0.05n c (equivalent to 5 x 10 19 cm -3 , which is 
similar to the experiment). The lengths of the density ramps at the rear vacuum plasma 
interface were /im, 100 /im and 200 /im. The centre of each of these ramps was at the 
same position so that the total amount of plasma the laser travels through is the same for 
all of the simulations. The laser pulse was linearly polarised with the laser's electric field in 
the y-direction, with tl = 500 fs and A = 1.053 /im. It was focused to a w = 2.4 /im at 
the top of the front density ramp to give an a ~ 15. 

The simulations assume an initially fully-ionised plasma and, in the case of an interaction 
with a relativistic laser pulse, the effects should be minimal as the front of the laser pulse will 
rapidly ionise any atoms it passes by. Care must be taken, when examining the simulation 
results when periodic boundary conditions have been used, to make sure the effect of any 
recirculation of particles is minimal, or to look at the result before recirculation occurs. The 
boundary conditions for these simulations were periodic in the transverse direction. 
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A. Ion acceleration and collimation mechanism 

After 1 ps into the simulation, the laser pulse starts to emerge from the rear of the plasma 
into vacuum. The accelerated electrons also move out into the vacuum region at the rear of 
the plasma. The ions do not immediately respond to the movement of the electrons due to 
their large mass and therefore a large space charge electric field in set up. The electron and 
ion densities can been seen in fig. [2] (a) and (b) at a time of 1 ps into the simulation, which 
shows the charge separation. The longitudinal electric field in the simulation, E S [ m is at its 
maximum value of « 0.7 TV/m at 1 ps and is shown in fig. [2] (d) (i). The simulation ion 
and electron data can be used to calculate the charge density and therefore the longitudinal 
electric field due to charge separation, E p , can be found. The charge separation contribution 
to the longitudinal electric field at 1 ps into the simulation is shown in fig. [2] (d) (ii) and 
has a maximum value of ~ 0.5 TV/m. Although the charge separation makes up the largest 
contribution to the longitudinal electric field, there is a discrepancy so there must be an 
additional electric field generation mechanism. 

It has previously been suggested that magnetic field effects will be important in similar 
regimes [ZU [22], [23J, [24] . The electron current generates a magnetic field around it and, as it 
moves into the vacuum, the magnetic field follows. Magnetic flux is conserved as the plasma 
expands into the vacuum and this can be seen in fig. [5] (c). This time- varying magnetic 
field contribution to the longitudinal electric field, E^ can be calculated by studying the 
temporal evolution of the vector potential due to the quasistatic magnetic field alone using 
V x B = — V 2 A. The induced electric field can then be calculated from Eg = — ^ and this 
contribution at a time of 1 ps into the simulation is shown in fig. [2] (d) (iii). The peak value 
of Eg at this time is ~ 0.2 TV/m and this makes up the discrepancy in the longitudinal 
electric field seen in the simulation, E sim . Furthermore, the contribution E p is larger than 
Eg at all times in the simulation. 

After 1.5 ps, a strong radial (^-direction) electric field is also observed at the exit of the 
channel, as can be seen in fig. [2] (e). This radial electric field acts to focus ions. The effect 
of the focusing force on the ions can be seen in fig. [2] (f), where the expanding ion front 
contains structure in the central region. By analysis of the charge density, it is possible to 
determine that the radial focusing electric field forms because of charge separation. The 
quasistatic magnetic field pinches the electrons, which therefore produces the collimating 
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FIG. 2: (a) The electron density, (b) the ion density and (c) the azimuthal magnetic field at a 
time of 1 ps into the simulation, (d) The longitudinal electric field components at a time of 1 ps 
into the simulation: (i) the total longitudinal electric field seen in the simulation, E s i mi (ii) the 
charge separation contribution, E p and (iii) the time varying magnetic field contribution, Eg. (e) 
The radial electric field at 1.5 ps, (f) the ion collimation at 2.5 ps and (g) shows the He 2+ ion p x 
against p y momentum at 2.5 ps. 
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FIG. 3: He 2+ ionp x against x-direction at times of 1.5 ps, 2.5 ps and 3.5 ps. The initial boundaries 
of the rear density ramp are shown by the dashed lines. 
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electric field for the ions. 

The ion momentum, p x against p y is shown in fig. [2] (g) and clearly shows the beam of ions 
accelerated in the x-direction and the radially accelerated ion populations. The transverse 
ion momentum, p yj is due to the radial ponderomotive Coulomb explosion and begins early 
in the simulation. Only once the electrons leave the target and the back surface sheath field 
is set up, the beam of ions is seen in the longitudinal (x) direction. The He 2+ ion beam has 
a maximum energy of about 45 MeV and a divergence of 7° half angle at 25 MeV at the end 
of the simulation. 

The ion acceleration region is illustrated in fig. [3j which shows the longitudinal momen- 
tum, p x , against the x-direction for times into the simulation of 1.5 ps, 2.5 ps and 3.5 ps. 
The dashed lines indicate the initial rear density ramp boundary. From this it can be seen 
that all of the high energy ions accelerated in the forward direction originate from the back 
of the density ramp near the plasma vacuum interface. It is interesting to note that the 
highest energy ions are accelerated from a significant depth into the density ramp. This 
suggests that the Debye length in this region is ~ 10 /im. 

The azimuthal magnetic field will exert a magnetic pressure, Pg = £> 2 /2/xq, on the 
electrons, preventing the return current, which will allow the electric field to persist for 
longer. Also, this will act to enhance the accelerating field for the ions by preventing cold 
electrons from entering the electric sheath region and allowing the ions to feel the accelerating 
fields for longer [211 1221 1231 El] . Holding back the background electrons will keep the electron 
temperature in the sheath region high, and the electron density low and thereby extending 
the Debye length in the region. It is difficult to quantify the overall influence of the magnetic 
pressure on the ion acceleration process from the simulations. But it can be deduced that 
the magnetic field presence is not a requirement for ion acceleration, since ID simulations 
show similar ion acceleration for a purely space charge generated electric field. 

B. Effect of plasma density 

A number of simulations were performed to investigate the influence of the plasma density 
on the ion acceleration. The density scan reveals the same dependance as that observed 
experimentally. At the lowest density (0.01n c ), the ion acceleration is almost purely in 
the transverse direction. However at 0.05n c the longitudinal acceleration is already more 
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density (n ) 

FIG. 4: Maximum He 2+ ion energies in both the longitudinal and transverse directions at 1.9 ps 
into the simulations, the maximum electric field (not at the same times in the simulations) and 
the total longitudinal energy in the electrons (> 10 MeV) for each of the densities. The lines are 
a smoothed curve fit to the data. 

effective than the radial acceleration. Fig. [4] shows the maximum He 2+ ion energies reached 
in the transverse and longitudinal direction at 1.9 ps into the simulation as a function of 
simulated plasma density. In the transverse direction, the ions are accelerated beyond the 
ponderomotive potential, U p = m e c(y/l + (a 2 ) — 1) « 5 MeV for peak a = 15, for all 
densities simulated. This can be attributed to self- focusing of the laser (enhancing a ) p3] 
and shock acceleration, as previously reported by Wei et al. [15]. 

The longitudinal acceleration increases with density until it peaks at 0.75n c , before a 
slight reduction as the density is increased to 0.9n c . As the density increases, there are 
more electrons available for acceleration and therefore density of exiting electrons on the 
rear surface increases. In figure [4] this is illustrated by the maximum electric field seen in 
the simulation and the total longitudinal energy in the fast (> 10 MeV) electrons, which 
can be seen to increase with density. This leads to an increase in the electric field and is 
able to accelerate ions to higher energies. 
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C. Effect of plasma ramp length 



To investigate the effect of the plasma ramp length, simulations with rear density ramp 
lengths of /xm, 100 /xm and 200 /im were performed. Due to the different positions of 
the plasma vacuum interface for the different ramp lengths, the ions start to accelerate at 
different times into the simulation. The simulations are compared at the time when the ion 
front has moved 100 /xm from the initial rear plasma vacuum interface. For the /im ramp 
this is at 3.2 ps, for the 100 /xm ramp this is at 3.9 ps and for the 200 /im ramp this is at 
4.4 ps. Fig. [5] shows the longitudinal ion spectra (left) and half-width angular divergence 
(right) of the ion beam as a function of energy at a time when the ion front has expanded 
by 100 /xm. The maximum electric sheath fields seen in the simulations are; 2.4 TVm -1 at 
a time of 1.7 ps for the /im ramp, 1.3 TVm -1 at a time of 2.2 ps for the 100 /im ramp and 
1.2 TVm" 1 at a time of 2.4 ps for the 200 /im ramp. As expected, the shortest density ramp 
allows the acceleration of the highest energy ions. However, it does not produce the most 
collimated ion beam, with the 100 /im and 200 /xm ramp lengths showing a significantly 
better collimation in the high energy tail than the /xm ramp. 

To investigate the improved collimation for longer ramp lengths, we need to consider the 
azimuthal magnetic field. Fig. [6] shows the azimuthal magnetic fields at a time of 2.2 ps after 
the start of the simulation for each ramp length. The longer ramps sustain a longer and 
stronger azimuthal magnetic field, which will increase the pinching of electrons and provide 
a stronger focusing force for the ions. The maximum magnetic fields seen in each of the 
simulations are 88 MG for the /xm ramp at 2.1 ps, 102 MG for the 100 /xm ramp at 2.0 ps 
and 85 MG for the 200 /im ramp at 2.3 ps. Even though the maximum field is slightly lower 
for the 200 /xm ramp, the energy stored in the magnetic field is greater as it is sustained 
over a longer length. 

IV. DISCUSSION 

A. Comparison with sheath acceleration models 

Consider a simple plasma expansion into a vacuum model of ion acceleration by an 
electrostatic sheath [25l [26] , which gives an expression for the maximum ion energy, e max ~ 
2Z/c B T e [in (r + a/t 2 + l)] , where T e is the fast electron temperature, r = u p it/\/2eE, t is 
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FIG. 5: The simulated He 2+ ion spectra (left) and the half-width divergence angles (right) as a 
function of energy for the different ramp lengths at a time when the ion front has expanded by 
100 fim. 
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FIG. 6: The azimuthal magnetic field at a time of 2.2 ps into the simulations for each of the ramp 
lengths. The initial rear density ramp positions are shown for each simulation. The simulation box 
is 598 x 168 /im. 



the time that the electrons are hot and e# denotes the numerical constant 2.71828 . . .. The 
model implies that the maximum ion acceleration achievable would improve with increasing 
the fast electron density (uj pi oc yjnl), temperature, T e or acceleration time, t. Taking a 
plasma density of 4 x 10 19 cm -3 with a measured electron temperature of 7.4 MeV and a 
laser pulse length of 1 ps, the model gives a maximum He 2+ energy of 117 MeV. However, the 
model assumes a sharp plasma vacuum boundary, which was not the case for the presented 
experiment and a density ramp is expected to reduce the maximum ion energy. 
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B. Laser parameter effects 



As the laser energy increases, it has been observed experimentally that the high energy 
electron spectra becomes hotter and contains more electrons [27J. The plasma expansion 
model therefore suggests that increasing the laser energy would be favorable to the ion ac- 
celeration. However, as the intensity is increased, the laser can break up into a number of 
filaments [28]. Due to the filament ation, there is a lower peak intensity but electrons can 
still be accelerated inside the individual plasma channels. This leads to the electron beam 
exiting the plasma in a number of lower temperature beams from a larger area. Due to com- 
putational constraints, the size of the simulation box used for the investigations performed 
here are not large enough to simulate the entire plasma length, which is presumably why 
filamentation of the laser beam did not occur. However, in similar parameter, but longer 
plasma, filaments are often observed. 

Although the effect of varying the laser pulse duration, r^, has not been studied here, 
assuming that the laser pulse length is directly related to the acceleration time, the plasma 
expansion model suggests that a longer pulse length would lead to higher energies. However, 
lengthening the laser pulse would not necessarily improve the acceleration since for fixed laser 
energy the intensity would decrease and give less efficient electron acceleration. Studies seem 
to imply that T e oc ao oc ^/rE [27] , whereas lengthening the pulse may lead to an acceleration 
time increase oc t. Therefore, there is likely to be an optimum pulse length for a particular 
laser energy, to maximize the ion energies and this has yet to be investigated. 

C. Plasma parameter effects 

Both the experiment and simulations show that the longitudinal ion acceleration improves 
with increasing ion density (to a certain point). It is expected that more electrons would be 
accelerated from a higher density plasma, leading to stronger accelerating sheath fields and 
this is shown in figure [4] as the total number of electrons accelerated. On the other hand, 
the electron temperatures have been observed experimentally to drop at higher plasma 
density [27] . which is possibly due to filamentation of the laser pulse. It is probable that the 
electron acceleration is overestimated in these 2D simulations as effects such as self-focusing 
are reduced and that is expected to lead to less filamentation. A reduction in the rate of the 
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acceleration improvement is seen in the maximum ion energy starting at 0.3n c , which could 
be related to the laser propagation. The group velocity of the laser, v g = cy/l — {nj (7) n c ), 
is slower in higher density plasma. This leads to the electron current outrunning the laser 
pulse. Not only does this reduce the energy the electrons can gain, but also the electrons 
arrive at the rear surface before the laser fields. Without the laser fields, the assistance of 
the ponderomotive force in removing the background electrons is not present and the sheath 
field would be reduced. This is a possible explanation for the reduction in maximum ion 
energy as the density approaches the critical density. 

For solid density targets, increased rear-side density scale-length, L, has been shown to 
dramatically reduce the acceleration, since the maximum electric field is given by E max ~ 
T e /{eL) [20, 29J. The Debye length, A^, for solid target plasma densities is on the order of 
nanometers. The scale-length due to the initial gas density ramp in both the experiment 
(~ 250 /im) and the simulations (71 /xm) is larger than the Debye length, which is estimated 
to be on the order of micrometers. Therefore, by the previous reasoning, the relatively long 
scale-length for the underdense plasma may be expected to reduce the maximum electric 
field by greater than an order of magnitude. On the other hand, the magnetic pressure 
on the background electrons, which are impeded from entering the sheath field region, will 
have the effect of keeping the Debye length in this region relatively long. Another interesting 
effect associated with the rear plasma ramp length is the improved collimation of the ion 
beam with increasing ramp length as seen in the simulations, which is again a consequence 
of the magnetic field generation. 



V. SUMMARY 



A collimated beam of ions was measured from an experiment performed using a high- 
intensity laser interacting with a helium gas jet target. Simulations have shown that the 
ion acceleration mechanism was target normal sheath acceleration with contributions to the 
electric field from both charge separation and the time- varying magnetic field. The influence 
of the time- varying magnetic field on the interaction and ion acceleration is expected to be of 
greater importance for ultra-short (t& ~ fs) petawatt class laser interactions. The maximum 
ion energy is expected to increase with plasma density due to the larger maximum electric 
field generated by the greater electron current. Indeed, the simulations suggest that very 
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high energy ions could be accelerated if sub-critical, but close to critical densities could be 
accessed. 
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